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The role of laser frequency chirps in the laser wakefield accelerator is examined. We 
show that in the linear regime, the evolution of the laser pulse length is affected by the 
frequency chirp, and that positive (negative) chirp compresses (stretches) the laser pulse, 
thereby increasing (decreasing) the peak vector potential and wakefield amplitude. In the 
blowout regime, the frequency chirp can be used to fine tune the localized etching rates at 
the front of the laser. In our simulations, chirped laser pulses can lead to 15% higher self- 
trapped electrons, and 10% higher peak energies as compare to the transform-limited pulse. 
Chirps may be used to control the phase velocity of the wake, and to relax the self-guiding 
conditions at the front of the laser. Our predictions are confirmed by multi-dimensional 
particle-in-cell simulations with OSIRIS. 



I. INTRODUCTION 



One of the main goals of plasma based accelerators is to deliver multi-GeV electrons [IHl] in 
distances that can be orders of magnitude shorter than with standard acceleration techniques. 
Recent particle-in-cell simulations in Lorentz-boosted frames [4, 5j predict electron bunch energies 
beyond 10 GeV in meter-scale plasmas using next generation 10 Petawatt laser systems. In fact, 
accelerating wakefields of nearly 50 GeV/m have been observed experimentally in plasmas [1], 
which are almost 1000 times higher than the fields observed in conventional accelerators. In 
laser or plasma wakefield acceleration (LWFA/PWFA), a short laser pulse or ultra relativistic 
electron beam propagates in an underdense plasma, and excites plasma waves [H [7] that can trap 
and accelerate electrons to ultra relativistic energies. This paper examines the role of frequency 
chirped laser in the LWFA, where the laser pushes plasma electrons away from the propagation 
axis through the ponderomotive force. This creates a positive space charge, as the ions remain 
essentially immobile in the time scales associated with the plasma period. Depending upon the 
laser and plasma parameters, linear or nonlinear plasma waves are excited. If the laser intensity 
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is sufficiently high, the radial ponderomotive force can lead to the cavitation of all the plasma 
electrons from the region where the laser propagates, creating a spherical plasma wave (bubble or 
blowout regime [SI [U] ) . 

The potential of the blowout regime for several applications has been confirmed in numer- 
ous experiments. The large accelerating fields associated with the blowout, which can trap and 
accelerate plasma electrons (self-injection), lead to the generation of quasi- monoenergetic multi- 
GeV electrons [2t [3l [T0Vll2j . In addition, the linear transverse focusing forces associated with the 
bubble [9] are ideal for the generation of X-ray radiation, as the accelerated electron beams per- 
form betatron oscillations in the ion channel |13H16| . The beams obtained from LWFA also have 
potential to drive a free electron laser [T7] after solving their present beam quality issues. For 
these applications it is crucial to control and manipulate the injection process that determines the 
charge, energy, energy spread, and strength parameter for X-ray radiation [18j. Although several 
techniques have been proposed to this end, including the use of short plasma down-ramps \19\ 120], 
the use of transverse external magnetic fields [2lj , and through the beating structures associated 
with counter or cross propagating lasers [21 [22[ [23] , and through ionization mechanisms [24J , this 
paper explores the possibility of using chirped lasers to control self-injection. 

Previous investigations on the role of the laser envelope asymmetries in wakefield excitation 
have already shown that a sharp laser intensity rise can drive stronger wakefields |25H27j . In 
addition, theoretical and simulation work |28l [29] on the impact of frequency chirps on the long 
pulse instabilities have shown that the growth of Raman forward scattering like instabilities can 
be controlled by acting upon the laser frequency chirps. The role of frequency chirps on the laser 
intensity profile has been explored experimentally [30] in the self-modulated LWFA, showing a 
significant enhancement in the total charge for sharp rising asymmetric pulses with positive chirps, 
where the frequency is lower at the front of the laser pulse. In this case, the ponderomotive force is 
also stronger at the front of the laser, which leads to stronger wakefields. Supported by analytical 
results, these experiments then emphasized the importance of the asymmetry of the laser pulse on 
the wake excitation and particle acceleration. 

In this paper we investigate, through numerical simulations with 2D and 3D particle-in-cell 
(PIC) simulations in OSIRIS [31j, the role of the frequency chirps in the blowout regime to show 
that the chirp can be used to adjust the self-injection rates, charge, and the output energy of LWFA. 
For the same plasma density, even though the peak energy of the accelerated electrons are in the 
same level for the chirped and un-chirped laser pulses, the total injected charge can be increased 
by 15%, using a chirped pulse laser for these parameters. Positive chirps increase (decrease) the 
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self- injected charge (maximum energy) by up to 10% in comparison to negative chirps in state-of- 
the-art conditions. Moreover, frequency chirps also change the laser group velocity in the plasma. 
These results agree with Ref. [30], although they are due to a different physical mechanism. 

In Section |IIJ an analytical model for the longitudinal bunching of a laser pulse with a frequency 
chirp is developed in the linear regime. Our model shows that laser with positive (negative) 
frequency chirp will compress (stretch) throughout the propagation. Accordingly, in the linear 
regime positive (negative) chirp leads to higher (lower) peak laser intensities. Good agreement 



between the analytical model and the simulations is found in the linear regime. In Section III 
the evolution of linear wake is investigated for different frequency chirped driving laser pulses. 
Although initially the wakefield amplitude is nearly independent of the sign of the frequency chirp, 
at later stages of the laser propagation, the wakefield amplitude increases (decreases) for positively 
(negatively) chirped pulse. In the blowout regime, which can not be examined by the linear model, 
the peak intensity increases for positive as well as negative chirps; however, the positively chirped 
pulse evolves faster than the negatively chirped pulse. In Section [TV] the effect of frequency chirp 



on self-injection in the blowout regime is analyzed. Finally, the conclusions are stated in Section [V] 

II. LONGITUDINAL BUNCHING 

The dominant contributions to the plasma refractive index, associated with the propagation of 
short laser pulses, are due to the ponderomotive and effective mass nonlinearities. The nonlinear 
coupling between the laser pulse and the plasma dynamics can lead to several processes, such as 
laser self-steepening [27] , self-compression [32] , self- focusing [33] , and self- modulation [TJ [33] . For 
our work, however, the most relevant mechanisms are associated with the self-compression, by 
which the laser pulse length changes during its propagation in the plasma. 

In order to examine the self-compression of a chirped laser pulse, we will use similar physical 
arguments to those presented in [33] . Our analysis is valid as long as the envelope approximation 
can be employed, i.e., as long as koLQ >> 1, where ko is the central laser wave number, and Lq 
is the laser pulse length. This assumption is well verified in the LWFA, where the refractive index 
for a linearly polarized laser pulse, in the weakly relativistic regime, becomes [33| . 



V 



(1) 



2ujq [ "^0 2 ujQ 

and where oOp = \J noe^/ (meo) is the plasma frequency, no is the background plasma density, 
Wo is the laser frequency, a = eA/(meC^) is the normalized vector potential, me is the electron 
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mass, c is the speed of light, and where 5n = n — uq and Su} = 00 — ujq are perturbations to the 
plasma density and laser central frequency ujq. In addition, (a^) is the average of over one laser 
period. According to the Eq. the laser group velocity of the laser Vg = {c'^/v^)/[l + iOpi^jo — 
l)/(2a;Q7o(7o + 1))] tSl El], where V(j, = cr]~^ is the phase velocity and 70 = y^(1 + (a^)/2), can 
be written as [33] 



l_^/l + ^_^_2 — 1 



2u}q \ ~^ no 4 ujQ j 



(2) 



In Eqs. Q and ([2]), the second term within the brace ({}) is related with the plasma density 
perturbations 5n, the next term represents the relativistic effects associated with the quiver motion 
of the electrons in the laser field, and the last term is due to the laser frequency modulations. The 
evolution of the laser pulse length is given by |33j 

Ldt c dC ^ ' 

where L is the laser length, t is time, and ^ = t — z/c is the distance in the co-moving frame. Using 
Eq. Q, Eq. ^ can be rewritten as 



dlogL _ ojI d 
dt ~ 2^8^ 



6n {a?) ^5uj 
no 4 Wo 



(4) 



To examine the relative contribution of these terms, we can consider, in the linear regime, the 
laser pulse profile (a^) ~ aosin^(c7r^/L). For this profile, the density modulation in the wake varies 
as (5n/no ~ 7raQsin(a;p^)/8 [7], and the maximum contribution from the first two terms on the 
right hand side of the Eq. Q is ~ ujpa^/2. In the regime where 5oj/uiq » LOpO^/A, the effect 
of ponderomotive and relativistic nonlinearities can be neglected with respect to the term for the 
spectrum change. In an underdense plasma (wp/wg << 1), for a gaussian chirped laser pulse with 
electric field E = ^oexp {—c^^^/L?) exp(— «a;o(l + /S^O ^"^^ chirp coefficient /3, the average wave 
number {k) = f^^kWdk/{f^^Wdk) [MlIM]; where W = f dsE{^ - s/2) ■ E*{^ + s/2) exp{iks) is 
the Wigner transform, can be written as (fc) = A;o(l + 2/3,^), with ko ~ ojq and 6uj/ujq ~ 5{k) /k^ = 
2/3^. Therefore the first and second term on the right hand side of Eq. Q can be neglected when 
2/3^ ~ 2/3L >> (jjpO^/L. In the remainder of this paper /3 >> ujpa'^/2 will be assumed. 

To evaluate further Eq. Q we note that in the absence of any nonlinear plasma 
effects, the frequency chirp coefficient of a gaussian chirped pulse varies as /? = 
[2c^ ln2/(L^CL!o)]-\/(-L/Lm)^ — 1) [37J, where Lj„ is the pulse length of the transform-limited laser 
pulse. As the laser pulse disperses with time, the rate of frequency variation inside the pulse 
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FIG. 1: Effect of the initial laser frequency chirp on longitudinal bunching: (a) variation of L (normalized 
by Lq) with time t (normalized by OJ^^) in the linear regime (aq — 0.05, uio/iOp = 20, Lq = 'Sc/cup), and (b) 
L variation with time in the nonlinear regime (aq = 2.0, LOo/ujp — 20, Lq ~ bc/up) for various /3o, ranging 
from — O.Olojp to O.OlWp. Solid lines in (a) show the theoretical prediction for the pulse length evolution 
with time for different chirp, (c), (d), and (e) plot the Wigner transform of the positively chirped laser 
(/3o = O.OSwp) used in the linear regime (a), and solid lines represent the (k) distribution within the pulse 
at u!pt ~ 0, 500 and 1000 respectively. 
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changes, and the frequency chirp coefficient varies as /3 = (3o{Lq / L) y {L^ — L^^) / {Lq — L^^) , 
where Lq and /3o are the initial pulse length and the frequency chirp coefficient respectively. For 
L, Lq » Lin, frequency chirp coefficient varies as /? = PqLq/L. Equation Q then yields: 

1/2 

(5) 



L = Ln 




1_ - i_2/3ot4— ^ + 



m 



where Lm/Lo = + (LqUJoPo)^ / {2c'^\n2)^]. For the typical parameters of interest for LWFA 

(Lo ~ Ap/2, loq = 20LOp; Xp = 27rc/ujp is the plasma wavelength), and fio/uip ~ 0(10"^), 
{Lin/L^f ~ 0.37 , in this case the expression for the pulse length [Eq. ([S])] simplifies to 
L = Lq{1 — 2(3QtiOp/ujQ), which predicts that the positive (negative) frequency chirp compresses 
(stretches) the laser pulse as it propagates through the plasma. This can be interpreted by inves- 
tigating the dynamics of different laser photons. For /3o > (/3o < 0), the photons located at the 
front move slower (faster) than the photons located at the back (front) of the laser. Thus, the 
distance between the laser photons decreases (increases) as the laser propagates, and laser length 
thus decreases (increases). 

In order to further investigate the role of the frequency chirp in the laser pulse length evolution, 
ID OSIRIS PIC simulations j31j were performed. The simulation uses a moving window that 
travels at c, with length 30 c/up, and divided into 3000 cells, with 500 particles per cell. The 
length of the plasma is 1000 c/ojp, and the ions form an immobile neutralizing fluid background. 
In this section the role of the chirp is identified by keeping the laser length, oq, and ujo/ujp constant 



regardless of the amount of chirp used. In Section IV, studies are shown, where we change the 
pulse length and the vector potential of the laser consistently with its chirp. 

The prediction for longitudinal compression [Eq. ([5])] is in good agreement with the simulations, 
as illustrated in Fig. [ij Figure [T|^a) uses a laser with ao = 0.05, ujo/ijjp = 20 and Lq = 3 c/ujp, thus 
exciting linear plasma waves. The laser pulse compresses (stretches) linearly with time for positive 
(negative) frequency chirp. We note that in these conditions, and in comparison to the scenarios 
where /3o 7^ 0, the pulse length remains constant during the laser propagation for /3o = O.Owp 
(if dispersion effects are neglected). This thus indicates that the ponderomotive and relativistic 
nonlinearities cancel, confirming that the pulse compression is essentially determined by the ini- 
tial frequency chirp. Figures [T]^c), (d) and (e) show the Wigner transform j35l [36], and average 
wavenumber distribution {{k)) within a chirped laser pulse with /3o = 0.05 Up at Upt = 0, ojpt = 500 
and ojpt = 1000. The variation in average wavenumber is linear within the pulse until ujpt ~ 500, 
and varies as {k) = A;o(l + 2/3^), where /3 = /3o = 0.05a;p at Upt = 0, and /3 ^ (3q(L/Lq) = 0.057 ujp 
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at ojpt = 500. For time ojpt > 1000, (k) does not change hnearly within the laser pulse. The 
laser evolution can not be predicted by the theory discussed here for ujpt > 1000 since non-linear 
mechanism such as self-steepening [27j starts to be relevant. 



Figure [T|^b) shows the evolution of the laser length using oq = 2.0 for various frequency chirps 

varying from /3o = —0.05 cOp to 0.05 LOp. In this particular case, pulse compression is observed for 
the range of frequency chirps (3q > —0.05 cOp] however, the compression was relatively stronger 
for positively chirped pulses. At higher intensities (oq ~ 1), in conditions where our analytical 
model is not valid, the two nonlinearities (ponderomotive and relativistic mass) do not cancel each 
other entirely, resulting into the net compression of the pulse [27], and including the contribution 
of frequency chirp will produce the effects as shown in Fig. fltb). 




In this section we have shown that the initial frequency chirp can compress or stretch the laser 
pulse depending upon the sign of /3o- In the next section we examine the effect of frequency chirp 
on the Wakefield excitation in the linear regime due to the longitudinal bunching. 



The plasma density modulations driven by a linearly polarized laser with normalized vector 
potential = ao x 'Z'{z,t) x cos[a;o(l -|- are given by [71 [38] 



where, oq is the initial peak normalized vector potential, Z{z, t) is the longitudinal profile of a^, his 
the density modulation in the wake, and uq is the initial homogeneous plasma density. Considering, 





III. WAKEFIELD EXCITATION IN LINEAR REGIME 




(6) 



the laser electric field in ID is given hy E = xEq ex.p[—^'^ / {2L{t)'^)] cos(a;^), then Z is given as 





where L(t) is given by Eq. Further assuming that /3L/c << 1, the ID solution for Eq. ^ 

is OEl] 
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FIG. 2: Evolution of peak plasma density (n,nax) as a function of propagation distance for different chirps 
and using gq — 0.05, ujq/lUp — 20, Lq — 3c/ujp. (a) shows the simulation values of TT-max at different times 
for the chirp coefficients /3 = 0.0 ojp by 'o', — O.OSwp by '-f', and O.OdUp by and (b) compares the 
theory(dashed line) and simulation(solid line) by plotting the rtmax normalized by ng with frequency chirp 
coefficient /3o (normalized by ujp) at early stage, neglecting the effect of pump evolution. Solid lines in (a) 
shows the theoretical predictions. 



In order to directly compare the theoretical predictions of the plasma density modulations with 
the simulations, Eq. ([s]) was solved using the same laser pulse profile as used in the simulations. 
In the simulations we have used a gaussian like 5th order symmetric polynomial profile defined as 

Z = 10/(z)3-15/(z)4 + 6/(z)^ (9) 

where, for —Lq < z < 0, f{z) = (Lq + z)/Lo, and for < z < Lq, f{z) = {Lq — z)/Lq. We used 
the simulation parameters from Section |ll| Figure [2]^b) shows the initial wakefield amplitude for 
different frequency chirps. For the typical values of /3Lo used in our simulations, no significant 
variation in the wake, at initial time, is observed with different /3s, e.g., for /3o = ±0.05 cjp, the 
wakefield amplitude changes by 1% of the wake amplitude driven by an un-chirped laser pulse. This 
is because, if we ignore the pulse length variation, the normalized peak laser vector potential is 
Oi/a^=o = 1/(1 — /3qLq/c^) ~ 1, being a^=o the peak laser vector potential for /3o = 0, no significant 
change in the wakefield with frequency chirp should be expected. Figure [2]^a) demonstrates, how- 
ever, that the frequency chirp (/3o > 0.05 ojp) can significantly change the amplitude of the plasma 
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wave during the laser propagation. Since the laser pulse compresses or stretches according to the 
sign of /3o (cf. Sec. |ll]), increasing or decreasing the laser vector potential, the amplitude of the 
plasma density modulation also changes accordingly. The evolution of the wakefield amplitude, 
plotted in Fig. ^h), shows that this amplitude increases (decreases) for positively (negatively) 
chirped pulse, as the laser peak vector potential increases (decreases). Inserting Eq. ^ in Eq. ([T]), 
the vector potential can be approximated as 

« alii + 2^otcol/4>xp[-f/Ll]/{l + M?, (10) 

which, on using in Eq. gives n w (1 + 2f3otUp/uQ)hQ, where uq is the initial wake amplitude. 
Thus, the density perturbation depends upon the chirp as 

n/no - 1 ~ 2/3oL^2/^2^_ (11) 

The theoretical estimates for the density modulations match with the simulation results as shown 
in Fig. [2j 

In a transform limited 20fs Ti:Sapphire laser pulse, the maximum chirp coefficient that can 
be introduced is around Pq ^ 1.0 x 10^^ sec~^ ~ O.Olcjp by stretching the pulse to ~ 28.3fs, for 
Wo = 20u)p. In such scenarios the wake amplitude changes by only 0.5% at Upt = 1000. Therefore, 
pulse compression, and hence the wakefield enhancement due to frequency chirp, play a significant 
role in LWFA for higher values of frequency chirp coefficients /3o > O.OSwp, which can be achieved 
by stretching a ~ 12 fs transform limited laser pulse to a chirped ~ 18 fs laser pulse. 

In the blowout regime, where the plasma dynamics is highly nonlinear, effects like self steepening 
and localized laser absorption also play a significant role in driving a nonlinear wake, and in self- 
injection. In such scenarios, frequency chirp may influence the laser etching rates, which further 
can affect the injection rates and beam characteristics. This will be analyzed in the following 
section. 

IV. BLOWOUT REGIME 

In order to investigate the role of laser frequency chirp in the blowout regime, a set of 3D PIC 
simulations were performed. For this purpose we consider a linearly polarized 350 mJ transform- 
limited laser with pulse duration 20 fs [full width at half maximum of the field(FWHM)], and 
6 ^m spot size, with central laser wavelength Aq = 800 nm. In the simulation, these parameters 
are translated into a transform- limited pulse with ag = 5.0, with central frequency ujq = 8 UJp, 
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FIG. 3: Time evolution of (a) peak energy, (b) total charge injected inside the first bubble, and (c) bubble 
radius in four scenarios (cases). Case (i): ['solid line' (black)] transform-limited pulse with ao = 5.0, 
luq = 8 LJp, Wq — 5 c/wp and L = 5 c/ujp. Case (ii): ['dash-dot line' (green)] same transform-limited pulse 
with flo — 5.0, loq = 11.3 Wp, Wq = 3.54 c/ujp and L — 3.54 c/cop. Case (iii) ['dashed line' (red)] positively 
chirped laser pulse with /3 ~ 0.0055ajp, ao — 4.23, loq = 11-3 Wp, Wq — 3.54 c/cup and L — 5 c/ujp. Case (iv) 
['dotted hue' (blue)] negatively chirped pulse with /3 = — 0.0055a;p, oq = 5.0, loq = 11.3 ujp, Wq = 3.54 c/ojp 
and L = 5 cjujp. Dark dots in Fig.[3](a) represents the time at which peak energy, as well as peak efficiency 
is achieved in LWFA. For cases (i), (iii), and (iv) the efficiencies are the same but greater than the efficiency 
for the case (ii) by ^ 15%. 

pulse length Lq = 5 c/wp, transverse spot size Wq = 5 c/wp, and with Gaussian transverse field 
profile given as exp(— r^/M^Q ), where r is the transverse coordinate, for a homogeneous plasma with 
plasma-density no = 1-75 x 10^^ cm~^. As the pulse is stretched, the maximum chirp coefficient 
f^raax = c^ln^/L?^ Can be obtained at L = 1.414 Lj.„ [37]. This corresponds to a chirp coefficient 
of Po = 0.0055 LOp, the peak normalized vector potential is reduced to oq = 4.23. Adding this 
frequency chirp stretches the pulse to ~ 30 fs. To keep the pulse length and plasma wavelength 
ratio equivalent to the transform-limited laser case, i.e., L = 5 c/ujp, electron density is lowered to 
riQ = 8.77 X lO"*^*^ cm~^, which translates into ojq = 11.3 ujp, and Wq = 3.54 c/ujp for a chirped pulse. 
For a negatively chirped pulse with /3 = —0.0055 ujp the rest of the laser simulation parameters 
are equivalent to the positively chirped pulse (/3 = 0.0055 ujp) case. These chirps can be routinely 
introduced/controlled in the experiments with lasers [40J. In order to compare directly the effect 
of chirp on the LWFA a simulation using the transform limited pulse propagating in a plasma with 
density no = 8.77 x lO^^cm^^ (same as used with the chirped pulses) was also performed, using 
simulation parameters as ao = 5, Lq = 3.54 c/cop, Wq = 3.54 c/up, and ujq = 11-3 ujp. The laser 
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is initialized in a simulation window that moves with the speed of the light, and with dimensions 
30c/a;p x SGc/wp x 36c/(x)p, divided into 1800 x 180 x 180 cells. The 3D simulations used 2 particles 
per cell. 

Figure [3] shows the variation in electron beam energy and total charge in the first bucket due to 
the introduction of frequency chirp in a transform-limited pulse in the 3D simulations. We discuss 
here results for the four cases, (i) transform- limited pulse with uq = 5, ojq = 8 Up, Lq = 5 c/ujp, 
and Wq = 5 c/ujp; (ii) same transform limited pulse (but with plasma density same as used with 
the chirped pulse) with = 5, ujq = 11.3 Wp, Lq = 3.54 c/cop, and Wq = 3.54 c/cop] (iii) positively 
chirped pulse with f3 = O.OObbujp, ao = 4.23, ujq = 11.3 cOp, Lq = 5 c/ojp, and Wq = 3.54 c/up] 
and (iv) negatively chirped pulse with /3 = — 0.0055u;p, uq = 4.23, ujq = 11.3 ujp, Lq = 5 c/up, and 
Wq = 3.54 c/iOp. Two scenarios are identified. First, when the ratio of pulse duration to plasma 
wavelength is kept constant by changing the plasma density [Cases (i), (iii) and (iv)]; and second, 
when plasma density is kept same [Cases (ii), (iii) and (iv)]. In the first scenario, for the chirped 
pulse [(iii) and (iv)] the peak energy of the accelerated electrons reaches 1.5 times the peak energy 
that can be reached using the transform-limited pulse [(i)], however on longer time scales [Fig. [3] 
(a)]. The higher peak energy for chirped pulse is obtained at the expense of lower total charge 
[Fig. |3] (b)] as compare to the transform-limited (un-chirped) laser pulse. The main reason behind 
these differences is the higher ojo/ojp, i.e., lower plasma density in the case of the chirped pulse. 
Since the laser group velocity (vg), and hence the wake phase velocity v^, which play a key role 
in the self-trapping mechanisms in the LWFA |2H I24j . is lower for larger densities, the trapping 
thresholds [7] in case (i) is relaxed as compare to the cases (ii), (iii) and (iv). 

In the case of transform-limited laser pulse [case (i)], the peak energy (200 MeV) is achieved 
after ~ 200//m. At ~ 400/im, the laser becomes pump-depleted, and the transition to the plasma 
Wakefield accelerator (PWFA) was observed. Since the blowout radius is measured at the same 
longitudinal position in Fig. [3](c), it drops to zero when the laser pump-depletes. Simulations show 
that self-injection still occurs when the wake is driven by the laser self-injected electrons, which 
increases the total charge of the beam at ~ 500^m. Around ~ 700;um, the self-injected electrons 
can no longer sustain the wake, and are lost to the background plasma. Total acceleration length 
for LWFA in this case is ~ 220/im, which is three times lower as compare to the acceleration 
lengths observed for the chirped pulse. The dark dots in Fig. [3] (a) represent the points where 
the peak energy is reached. The efficiency is similar for all the three cases. For the chirped pulse 
the self-injection and acceleration process is relatively gradual as compare to the un-chirped pulse. 
The amount of injected charge for the chirped pulse can reach to 60% of the total charge obtained 
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with the transform-limited pulse [case (i)]. 

In the second scenario (keeping plasma density constant), the peak energy for the transform- 
limited pulse [case (ii)] is ~ 15% lower than in the case of negatively chirped pulse [ Fig. |3]^a)], 
and with total peak charge ~ 16% lesser than the charge obtained with the chirped pulse [Fig. [s] 
(b)]. Hence, using a chirped pulse in place of a transform limited pulse provides better efficiency, 
as well as in ~ 25% shorter acceleration lengths. 

Within the chirped pulses, the positively chirped pulse provides higher (lower) charge (peak 
energy) as compared to the negatively chirped pulse [Fig. [sj. Simulations suggested that the higher 
final electron energies in the case of negatively chirped pulse is due to the fact that negatively 
chirped laser propagates with higher group velocities, followed by the bubble with higher phase 
velocity. This increases the threshold for the initial 7 of the electrons which can be trapped, 
reducing the number of electrons, as well as electrons that can be accelerated to higher velocities. 
Beam loading effects f^lj also play a significant role in reducing the wake field for the positively 
chirped pulse due to the higher charge self-injection, where as for negatively chirped pulse the 
beam loading effects are not significant due to relatively lower charge injection as compare to 
the positively chirped pulse. These simulations showed that by adjusting the initial chirp of the 
laser pulse, the number of self-injected electrons, and self-injected electron beam energy can be 
controlled. For the laser parameters discussed here, the difference in total charge for the two 
chirped cases [(iii) and (iv)] can reach up to 15%, and in peak energy 10% of the difference can be 
achieved. These results highlights some of the advantages of using chirped pulsed laser in LWFA. 



According to the Section III, the wake field amplitude is not significantly changed for the 
frequency chirp discussed here. We believe that it is the combined effect of longitudinal bunching 
as well as localized etching of the laser which leads to differences in the LWFA by positively and 
negatively chirped pulse. 

Considering the effect of localized etching [35], the laser group velocity Vg is given by 



Vn — Vetch, where Vg = c[l — a;p/(2a;(^etch) )] is the laser linear group velocity, 
2 



9 

Vetch = cu;'p/uj{^etch)''^ IS the etching velocity [42], and ^etch is the position of the localized 
etching within the laser in ^ coordinates. A positively chirped pulse, with a red shifted front, 
etches faster, and may thus propagate with lower group velocity, as compared to a negatively 
chirped laser pulse. Furthermore, the phase velocity of the plasma wave (wake), propagating 
at the back of a laser, is roughly equal to the group velocity of the laser; thus, for a positively 
chirped laser pulse the wake phase velocity is lower as compared to a bubble propagating behind 
a negatively chirped laser pulse. Electrons with sufficiently high velocity in the forward direction 
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(this velocity roughly equals the phase velocity of the bubble at the back) have higher probability 
of getting trapped inside the bubble [39], which means that the positively chirped laser pulse 
lowers the threshold for electron trapping by slowing down the wake. On the other hand, a 
negatively chirped laser pulse, will have a higher group velocity, thus increasing the threshold for 
trapping. 

To illustrate clearly the effect of positive and negative frequency chirp on localized etching and 
on the laser self-guiding [39], in Fig.|4| we show the results from 3D PIC simulations for two linearly 
polarized laser pulses with ao = 5.0, with central frequency ujq = 20a;p, pulse length Lq = 3c/ Up, 
and transverse spot size Wq = A.5c/ujp, and with /3 = ztO.OScjp. We have exaggerated the chirp 
such that effects attributed to the chirp are clearly visible. Since for /3 > the laser intensity is 
higher at the front than when /3 < 0, a sharper density spike is formed at the front of the laser 
for ;g > [Fig. Qa), (d)], which leads to a stronger localized pump- depletion [Fig. |4];b)-(f)]. For 
the negatively chirped pulse considered in the simulations, the density spikes are not sharp, so the 
localized etching is either absent, or very weak, as compared to the positively chirped pulse. Thus, 
positively chirped lasers may also relax the self-guiding conditions for stable laser propagation in 
the blowout regime. These results indicate that the variation in laser group velocity and localized 
pump depletion, due to the initial laser frequency chirp, may impact the rates at which self-injection 
occurs, and number of self-trapped particles. For the parameters considered in the Fig. |3] [cases 
(ii) and (iii)], we observe the effect of chirp sign (positive or negative) on localized laser etching; 
however, and for realistic chirps, the effects are not as strong as shown in Fig. [4j since the amount 
of frequency chirp used in the laser for Fig. |4] is 5 times higher as compare to the frequency chirp 
used in the laser for Fig. [3j 

V. CONCLUSIONS 

In conclusion, we have examined the effect of the frequency chirp on LWFA. In the linear 
regime, a positive (negative) frequency chirp compresses (stretches) laser pulse, resulting into 
enhanced (reduced) wakefield amplitude with time due to dispersive effects. In the blowout 
regime, and using the chirps that can be routinely introduced in the experiments, simulations 
show that negatively chirped lasers can provide higher peak energies to the self-injected electrons 
in comparison to un-chirped lasers. Moreover chirped laser pulses can also lead to higher number 
of self-injected electrons. In addition, the laser group velocity, and thus the wake phase velocity. 
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is lower (higher) when the laser frequency is lower (higher) at the front of the laser due to positive 
(negative) frequency linear chirp, which then influences the rate at which self-injection occurs and 
the number of self-injected electrons, providing an extra control over the self-injection process in 
LWFA. 
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FIG. 4: Effect of frequency chirp on localized laser etching: (a) and (d) line out of charge density on the 
axis at the front of the laser, (b) and (e) line out of laser field on the axis, and (c) and (f) laser field (Ex) 
in a; — z plane at time 448.43/a;p for /3 — 0.05 LOp{a, b, c) and /3 — —0.05 LOp{d, e, /) respectively. 
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